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ABSTRACT 

Aims. To study chemical composition of metal-rich absorbers at high redshifts in order to understand their nature and 
to determine sources of their metal enrichment. 

Methods. From six spectra of high-z QSOs, we select eleven metal-rich, Z > Zq, and optically-thin to the ionizing 
radiation, Ai'(Hl) < lO^'^ cm~^ , absorption systems ranging between z — 1.5 and z = 2.9 and revealing lines of different 
ions in subsequent ionization stages. Computations are performed using the Monte Carlo Inversion (MCI) procedure 
complemented with the adjustment of the spectral shape of the ionizing radiation. This procedure along with selection 
criteria for the absorption systems guarantee the accuracy of the ionization corrections and of the derived element 
abundances (C, N, O, Mg, Al, Si, Fe). 

Results. The majority of the systems (10 from 11) show abundance patterns which relate them to outflows from low 
and intermediate mass stars. One absorber is enriched prevalently by SNe II, however, a low percentage of such systems 
in our sample is conditioned by the selection criteria. All systems have sub-kpc linear sizes along the line-of-sight with 
many less than ~ 20 pc. In several systems, silicon is deficient, presumably due to the depletion onto dust grains in the 
envelopes of dust- forming stars and the subsequent gas-dust separation. At any value of [C/H], nitrogen can be either 
deficient, [N/C] < 0, or enhanced, [N/C] > 0, which supposes that the nitrogen enrichment occurs irregularly. In some 
cases, the lines of Mgll AA2796, 2803 appear to be shifted, probably as a result of an enhanced content of heavy isotopes 
^^Mg and ^^Mg in the absorbing gas relative to the solar isotopic composition. Seven absorbers are characterized by 
low mean ionization parameter U, log U < —2.3, among them only one system has a redshift z > 2 (zabs = 2.5745) 
whereas all others are found at z ~ 1.8. This statistics is not affected by any selection criteria and reflects the real rise 
in number of such systems at z < 2.0. Comparing the space number density of metal-rich absorbers with the comoving 
density of star- forming galaxies at z ~ 2, we estimate that the circumgalactic volume of each galaxy is populated by 
~ 10^ — 10* such absorbers with total mass < 1/lOOt/i of the stellar galactic mass. Possible effects of high metal content 
on the peak values of star-forming and AGN activities at z ~ 2 are discussed. 

Key words. Cosmology: observations - Line: formation - Line: proflles - Quasars: absorption lines 



1. Introduction quasar spectra. Individual element abundances, or abun- 
5^ dance patterns, derived from the metal-rich systems can be 
, The majority of metal-bearing absorbers detected in quasar ^^^^ reconstruct their chemical history and to estimate 
spectra demonstrate very low metalhcity of about 1/30 of ^^^^^-^^ contributions from SNe II, SNe la, and AGB-stars. 
the solar value, Zq, and below (Songaila & Cowie 1996; mi i • i • i- -i • /x^ 
Simcoe et al. 2004). Compared to the metal-poor absorbers, ^,^^1",™''°''', "^T""' 'P'""'' Tf^. distribution (Fan 
the systems with high metallicity ~ from solar to overso- f «/" ^OOl) and observations of high-redsduft galaxies 
lar - occur in quasar spectra quite seldom. A particular (Dickinson et al. 2003) show that the redshift z ~ 2 rep- 
interest to the metal-rich absorbers stems from the fact "^^^e^*^ ^ turning point m the global cosmic evolution. The 
that high metallicities are usually measured in the cen- '^tudy of the chemical enrichment at different z can lead 
tral parts of galaxies (Rich et al. 2007; Cohen et al. 2008; important clues concermng the factors which cause the 
ZoccaU et al. 2008; Davies et al. 2009) and are closely re- P^^^^ ^^^^^^ star-forming activities at z ~ 2. 
lated to the starburst activity and chemical evolution of To obtain chemical enrichment patterns of the ab- 
AGNs/QSOs (Hamann & Ferland 1999; D'Odorico et al. sorbing gas, accurate abundances of different elements 
2004; PoUetta et al. 2008; Silverman et al. 2009). At high are required. Abundance measurements in high-metallicity 
redshifts, the only way to gain an insight into the chemi- damped Ly-a absorbers (DLAs) and sub-DLAs where ion- 
cal composition of gas distributed over large cosmological ization corrections are supposed to be small are limited 
distances is through the analysis of absorption systems in either to low abundant species like Zn, Mn, and Cr, or to 
those which have transitions with low oscillator strengths 



* Based on observations obtained with the UVES at the VLT like Fe and Si since heavily saturated lines of other ele- 
Kueyen telescope (ESO, Paranal, Chile). ments prevent their accurate column density determination 
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(Prochaska et al. 2006; Pcroux et al. 2006; Meiring et al. 
2007; Meiring et al. 2008; Peroux et al. 2008). However, 
just Fe and Si (as a proxy for a-elements) which belong to 
the key species in the abundance pattern, can be severely 
depleted (Hou et al. 2001; Vladilo 2002; Vladilo & Pcroux 
2005; Quast et al. 2008), thus hampering significantly the 
identification of the enrichment agents. Note that metal- 
rich DLAs are very rare — up to now, only about a dozen 
is detected. 

On the other hand, in absorbers which are optically-thin 
to the ionizing radiation^ the ionization correction factors 
are large and, consequently, the resulting element abun- 
dances arc very sensitive to their values. In general, ioniza- 
tion corrections are determined by the ionization parameter 
U and by the spectral shape of the ionizing radiation. To 
increase the reliability of the measured abundances both 
a special selection of absorption systems and an adequate 
mathematical treatment of the inverse spectroscopy prob- 
lem are required. 

An approach employed in the present paper is based on 
the physical model of absorbing gas which assumes fluc- 
tuating gas density and velocity fields inside the absorber. 
This model is a generalization of a commonly used approx- 
imation of a plane-parallel gas slab of uniform density and 
a microturbulent treatment of the velocity field. The cor- 
responding computational procedure used to calculate the 
metal abundances includes also the adjustment of the spec- 
tral shape of the ionizing radiation (see Sect. 2). 

Since relative element ratios can differ from the solar 
pattern, it should be clarified on base of which species we 
conclude whether the gas is metal-rich or not. Mctallicity 
in stars is usually characterized by the iron content, that of 
strong absorption systems (LLSs and DLAs) — by the con- 
tent of zinc or sulfur which are not condensed significantly 
into dust. In optically-thin absorption systems, mctallic- 
ity is conventionally estimated on base of carbon. Thus, 
in the present metal-rich sample we included absorbers 
with carbon abundances near and above the solar value, 
(C/H)0 = 2.45 X 10-4 (Asplund et al. 2004). 

The selection of systems occurred in the following way. 
In the framework of our project to study the spectral en- 
ergy distribution (SED) of the ionizing radiation the quasar 
spectra were searched for the optically-thin absorption sys- 
tems revealing lines of different ions (Agafonova et al. 2005). 
The presence of at least one pair of ions of the same element 
(e.g. C II, C iv) was essential for the system to be included 
in the sample. These subsequent ions ensure that the mean 
ionization parameter can be estimated without any addi- 
tional assumptions. In total, about 50 absorbers identified 
in high-resolution spectra of ~ 20 QSOs with emission red- 
shifts 1.7 < z < 4.5 were selected and analyzed. The de- 
rived metallicities range from 1/1000 to several times solar 
values. A part of results was described in our publications 
dealing with the evolution of the background SED with 
redshift (Reimers et al. 2006; Agafonova et al. 2007) and 
the restoration of the SED of the outcoming quasar radi- 
ation on base of associated systems (Reimers et al. 2005; 
Levshakov et al. 2008, hereafter Paper I). These publica- 
tions included already several systems with high metal con- 
tents, however, enrichment scenarios were not elaborated. 
This is just what the present paper is focused on where 
we consider additional 7 metal-rich absorbers. With the 



^ Column density of neutral hydrogen A'^(Hl) < 3x10^^^ cm ^ . 



previously described systems, the sample of metal-rich ab- 
sorbers consists now of 11 systems with redshifts ranging 
from ^; = 1.5 to ^; = 2.9. 

2. Computational method 

Absorption systems are analyzed by means of the Monte 
Carlo Inversion (MCI) procedure described in Levshakov, 
Agafonova, & Kegel (2000, hereafter LAK), with further 
details in Levshakov et al. (2003). This procedure supposes 
that all lines observed in a metal system are formed in the 
gas slab with fluctuating density nii{x) and velocity v{x) 
fields (here x is the space coordinate along the line-of-sight 
within the absorber). Further assumptions are that within 
the absorber the metal abundances are constant and the 
gas is in thermal and ionization equilibrium. The proce- 
dure is applicable for optically-thin systems. The inputs 
are the observed line profiles and the ionization curves for 
each ion (hydrogen -I- metals) included in the analysis. The 
ionization curves (the dependence of the ion fraction on the 
ionization parameter U) are computed with the photoion- 
ization code CLOUDY version 07.02.01 (last described by 
Ferland et al. 1998) which in turn uses as the inputs an 
ionizing continuum and a set of element abundances. The 
fitting parameters (outputs) of the MCI procedure are: the 
mean ionization parameter Uq, the total hydrogen column 
density N^. the line-of-sight velocity dispersion, a^, and the 
density dispersion, ay, of the bulk material [y = nu{x)/no], 
and the contents Za of metals included in the objective 
function [Eqs.(29, 30) in LAK]. The fitting parameters are 
estimated by minimization of residuals between synthetic 
and observed line profiles. As in all Monte Carlo methods, 
integrals [here the intensity at a point within the line pro- 
file, see Eq.(12) in LAK] are calculated by modeling the 
distribution of the integrated fimction. This function is a 
convolution of the gas density and the local absorption coef- 
ficient (which depends on the velocity through the irregular 
Doppler shifts) and, hence, the integrated function is real- 
ized through the density and velocity distributions over the 
integration path. Both the distributions of nii{x) and v{x) 
are represented by their sampled values {rii} and {vi} at 
equally spaced points Xi along the line-of-sight, and optimal 
configurations of {n j and {vi} are estimated by the simu- 
lated annealing algorithm. With the estimated values of the 
fitting parameters and the distributions of nii{x) and v{x), 
we can calculate the column densities of all ions [Eqs.(24, 
25) in LAK] and the mean values for the kinetic tempera- 
ture, Tkin, the gas density, no, and the linear size, L, of the 
absorber along the line-of-sight (see below). 

If the metal abundanc;es Z,, obtained from the MCI pro- 
cedure differ from those used in CLOUDY to calculate the 
ionization curves, then these curves are re-calculated with 
the updated abundances and the MCI procedure is run 
again. The iterations stop after the concordance between 
the abundances used in CLOUDY and derived from MCI 
is reached. 

At this stage of calculations, the shape (SED) of the 
ionizing radiation (which comes into the MCI procedure 
through the ionization curves) is treated as an external 
parameter, i.e., it is keeping constant. In case when a 
first guessed SED cannot provide an acceptable solution, 
a complemented computational routine aimed at adjusting 
the spectral shape is applied. The corresponding algorithm 
based on the response surface methodology from the theory 
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Fig. 1. Hydrogen and metal absorption lines from the Zabs ~ 2.5745 system towards HE1347-2457 (solid-line histograms). The 
vertical axis is normalized intensities. The zero radial velocity is fixed at 2: = 2.5745. Synthetic profiles are plotted by the smooth 
curves. The wide profiles of the OviAAlOSl, 1037 lines show that the narrower line metal-rich system is embedded in a hot gas 
(seen also in Hi Ly-a and Ly-/3 lines as a difference between the observed and synthetic profiles). The insert in panel Cm A977 
demonstrates a non-zero residual intensity at the line center. The identified blend is marked in panel Si IV A1402. 



of experimental design is described in detail in Agafonova 
et al. (2005, 2007). It includes a parameterization of the 
SED by means of a set of factors and the estimation of 
their values in a way to ensure that all obtained physical 
parameters arc self-consistent and all observed profiles are 
described with a sufficient accuracy (i.e., ~ !)• 

It is known from both model calculations (Haardt & 
Madau 1996; Fardal et al. 1998; Madau & Haardt 2009) and 
from SEDs restored from observational data (Agafonova 
et al. 2005, 2007) that spectral shapes of the intergalac- 
tic ionizing background in the energy range E > I Ryd 
have complex forms with an emission bump at < 3 Ryd 
and a stepwise break at £^ > 3 Ryd. However, as noted 
in Agafonova et al. (2007; Sect. 3.1.2), at high metallic- 
ities the dependence of the ion fractions on the absolute 
abundances (especially on the contents of the most abun- 
dant elements such as carbon and oxygen) is comparable or 
even overrides the dependence on the spectral shape. For 
the systems considered in the present study this means that 
the inherent uncertainty in the values of the absolute ele- 
ment abundances makes it impossible to estimate the SED 
of the underlying ionizing continuum in full detail. That is 



why the spectral shape is represented and estimated in a 
simplified form including only 2 factors: a power law index 
a (j^~") in the range 1 < E < A Ryd, and a depth of a 
break at 4 Ryd. The power law index for the energy range 
beyond this break (E > 4 Ryd) was fixed as 1 — 0.8a. 

The mean linear size L of an absorbing cloud along the 
line of sight can be estimated as L = Nh/uq, where is 
the total hydrogen column density, and uq is the mean gas 
density of the absorber. The mean gas density is calculated 
from the mean ionization parameter Uq and the density 
dispersion ay (LAK): 



Here nph is the number density of photons (in units of 
cm~^ ) with energies above 1 Ryd 




where c, /i, z/c, and J,y arc the speed of light, the Planck con- 
stant, the frequency of the hydrogen Lyman edge, and the 
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Fig. 2. The ionizing spectrum corresponding to the ionization 
state observed in the ^abs = 2.5745 system (HE1347-2457) is one 
from the shaded area. For comparison the intergalactic ionizing 
spectrum at 2 = 2.5 modeled by Haardt & Madau (1996) is 
shown by the dotted line. The spectra are normalized so that 
Ji,{hv = 1 Ryd) = 1. The positions of ionization thresholds of 
different ions are indicated by tick marks. 

specific intensity (in ergs cm~^s~^ sr~^ Hz~^). J912 is the 
normalizing value and represents the intensity of the ioniz- 
ing radiation at 912 A (1 Ryd). As a reference, the value 
of Jgi2 for the intergalactic background at a corresponding 
redshift is taken from Haardt & Madau (1996). 

If a strong radiation source such as a quasar is located 
close (up to ~2 Mpc away) to the absorber, then the in- 
tensity J912 of the incident ionizing radiation exceeds sig- 
nificantly the mean intergalactic level (see, e.g., examples 
in Paper I). Thus, the gas density/linear size calculated 
on base of the intergalactic J912 represents in fact only a 
lower /upper limit, respectively. 

A few systems from the present study show signs of 
the incomplete coverage of the background light source. 
Computational treatment of such systems occurred as de- 
scribed in Sect. 2 of Paper I. The covering factor is assumed 
to be the same for the whole cloud, i.e. it does not depend 
on A (or, equivalently, on the radial velocity v). However, 
different covering factors are allowed for different ions. This 
follows from a model of the absorption arising in a gas cloud 
with fluctuating density. Namely, ions of higher ionization 
stages trace rarefied gas which can be quite extended (i.e., 
covers a greater part the background source) whereas low 
ionization ions originate in more dense and, hence, compact 
volumes (i.e., the covering factor is smaller). 

The laboratory wavelengths and oscillator strengths for 
all lines except Sim were taken from Morton (2003). For 
Sim Morton gives 1206.500 A (see comments on page 
190 in Morton 1991), however with this wavelength the 
Sim line comes out shifted by ~2.5 km s~^ relative to 
both Sill and Siiv lines. Since metal lines in absorption 
systems considered in the present paper are very narrow 
{FWHM < 15 km s~^ ) this shift becomes noticeable. On 
the other hand, other catalogs of atomic and ionic spectral 
lines (e.g., Bashkin & Stoner 1975; Kelly 1987) and catalogs 



of solar lines (e.g., Sandlin et al. 1986; Curdt et al. 2001) 
give for Si m the wavelength 1200.51 A. This value fits much 
better to what is observed in the absorption systems and it 
was therefore adopted for the present work. 

Solar abundances were taken from Asplund et al. (2004). 
Note that their solar abundance of nitrogen, (N/H)© = 
6 X 10^'', is ~1.4 times (0.15 dex) lower than that from 
Holweger (2001). 

Quasar spectra were obtained with the UVES/VLT in 
the framework of the ESO Large Program 'QSO Absorption 
Line Systems' (ID N0.I66.A-OIO6). Data reduction was per- 
formed by B. Aracil (Aracil et al. 2004). Except Q0329 385, 
all QSOs used in the present paper have been discovered in 
course of the Hamburg/ESO survey (Wisotzki et al. 1996; 
Reimers et al. 1996; Wisotzki et al. 2000). 



3. Analysis of individual systems 

3.1. Quasar HE1347-2457 
3.1.1. System at ^abs = 2.5745 

The system at z^hs = 2.5745 consists of the neutral hy- 
drogen Lyman series lines (Li-Lio) and many lines of dif- 
ferent metal ions (Fig. 1). Expected positions of the Fell 
AA2382.76, 2600.17 fines coincide with strong telluric ab- 
sorptions, Feiii A1122.52 falls in a noisy part of the spec- 
trum and cannot be extracted from the noise. The emis- 
sion redshift of the quasar HE1347-2457 is Zem= 2.578 
(based on Cii emission), so that the absorber is detached 
by only 300 km s^^ from the quasar. The profiles of the 
Ovi AA1031,1037 absorption lines differ from the profiles 
of other ionic transitions all of which exhibit similar line 
shapes. The observed intensities of low ionization (C 11, Si 11) 
and high ionization (C iv, Si iv) lines do not differ from each 
other significantly implying a rather low ionization param- 
eter Uo, log(C/o) < —1-7, for all types of the incident ion- 
izing spectra. At such Uq, a considerable amount of Ovi 
can hardly be produced. Thus, this absorption system orig- 
inates probably in a dense and relatively cold gas (seen in 
lines of Cii-Civ, Siii-Siiv, Niii, Mgii, and Aim) embed- 
ded in a hot and highly ionized medium seen in Ovi and 
Hi AA1215, 1025 lines. 

The absorption profiles of the C ill A977 and H i A1025 
lines^ have flat bottoms and reveal non-zero residual in- 
tensities at the line centers (0.02 and 0.01, respectively) 
pointing to the incomplete coverage of the background light 
source, i.e. the covering factor for C ill is C(C m)= 0.98 and 
for hydrogen C(Hl)= 0.99. Due to the detection of many 
lines of neutral hydrogen its column density can be esti- 
mated quite accurately: iV(Hi) = (9.5 ± 0.3) x lO" cm'^ . 
A clear (unblended) doublet of Civ AA1548, 1550 afiows us 
to estimate C(Civ)= 0.99 and iV(Civ) = (3.2 ±0.1) x lO^^ 
cm~^ . Because of the line blending, low S/N and/or the 
presence of a single unsaturated line only (like Sim A 1206 
and Nlll A989) accurate covering factors for other ions can- 
not be determined. Their limiting values are as follows: 
0.98 > C(Cii/Siii/Mgii) > 0.6, 0.98 > C(Siiii/Niii/Aliii) 
> 0.7, and 0.98 > C(Siiv) > 0.8. We note that covering 



^ The Hi A1215 line is black because of blending with Ly-Q 
from the highly ionized system. An appropriate A''(Hl) for this 
system is ~ 10^* cm~^ , i.e. its Ly-a is saturated but Ly-/3 is 
not. 
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# 


QSO 


2cm 




log 


Af(H), 


no, cm * 


L, pc 




Ref." 


1 


HE1347-2457 


2.578 


2.5745 


-2.3 


3.5E17 


>5E-3 


<25 


0.7 


1 


2 






1.7529 


-2.5 


(5 - 8)E18 


0.007 - 0.1 


12 - 500 


1.1 


1 


3 






1.5080 


-2.5 


(5 - 8)E18 


>7E-3 


<500 


1.1 


1 


4A 


HE0151-4326 


2.775 


2.4158 


-2 --1.5 


(3 - 7)E16 


>7E-4 


<30 


~1.0 


1 


AB 






2.4158 


-3.2 - -2.2 


(1 - 10)E18 


>(0.1 - l)E-2 


<300 


1.4 - 2.0 


1 


AC! 






2.4158 


_2 1.5 




^7F,-4 


\ -Lt-> 


±.\j 


\ 


5 






1.7315 


-2.3 


(5 - 6)E18 


> 5E-3 


<300 


0.9 


1 


6 


HEOOOl-2340 


2.26 


1.6514 


-2.6 


(5 - 6)E17 


>8E-3 


<20 


0.8 


1 


7 






1.5770 


-2.0 - -1.8 


(6 - 8)E16 


>lE-3 


<20 


~ 0.2 


1 


8 


HE0141-3932 


1.80 


1.7817 


-2.5 


(4 - 5)E17 


>7E-3 


<20 


1 


2 


9 


HE2347-4342 


2.88 


2.8980 


>-1.0 


~1E17 


>lE-4 


<300 


1 


3 


10 






1.7963 


-2.7 


~1E17 


> 2E-2 


<4 


0.6 


4 


11 


Q0329-385 


2.435 


2.3520 


~ -0.75 


~ 1E17 


>lE-4 


<300 


2-3 


3 



^Mean kinetic temperature in units of lO'' K. 

"References. — (1) present paper. (2) Reimers et al. 2005. (3) Levshakov et al. 2008. (4) Agafonova et al. 2007. 



factors may be slightly different for each of the two com- 
ponents comprising the absorption system under consider- 
ation. 

With such uncertainties in covering factors the shape 
of the ionizing spectrum cannot be restored uniquely. 
Nevertheless, some conclusions about the SED can be ob- 
tained. The incomplete covering of the background light 
source indicates that the system is located close to the 
quasar. Common first guess for the SED of the outcom- 
ing quasar radiation is a power law, i/^". However, pure 
power law spectra are not consistent with the observed 
line intensities. Namely, for all tried spectral indices a, 
—2.0 < a < —0.5, the mean ionization parameter Uq is 
low and the corresponding fraction of C IV is low as well, 
which makes the abundance of carbon, 

^ Civ Thi 



to be very high, 4 — 5 times solar value (and the abun- 
dance of magnesium even higher). With such a high metal 
content it is impossible to describe the observed profiles of 
hydrogen lines: in particular the central parts of the syn- 
thetic profiles of the A972.53 and A937.80 lines look too 
wavy following the pattern exhibited by metal ions. In or- 
der to decrease the metallicity, the ratio Thi/Tqiv should 
be lowered, and this occurs when the observed profiles of 
C ii/C IV and Si ii/Si iii/Si iv lines are described with higher 
values of Uq. As it was the case with the associated systems 
from Paper I, to fit to these conditions the ionizing spec- 
tra should have a break at 4 Ryd. The depth of the break 
depends on the adopted covering factors, but the intensity 
decrease of at least one order of magnitude at ii^ > 4 Ryd 
is required. On the other hand, a too deep drop of the in- 
tensity (the depth of the break > 1.3 dex) leads to the 
overproduction of Sim compared to the observed profile 
even at the lower limit of C(Siiii). A small difference be- 
tween the Civ/Cii and Siiv/Siii ratios makes preferable 
the power indices between 1 and 4 Ryd of a ~ 1.2 — 1.4. 

As already noted, at high metallicities (above solar) the 
relative element contents, especially those of the most abun- 
dant elements carbon and oxygen, affect quite strongly the 
calculated ion fractions (upon which the procedure of spec- 
tral shape estimation is based). Lines of different oxygen 



ions are not available in the Zabs = 2.5745 system, and 
the unknown oxygen content gives an additional degree 
of freedom in the procedure of the spectral shape estima- 
tion. These considerations restrict the acceptable spectral 
shapes as shown by the shaded area in Fig. 2. The values 
of the physical parameters obtained with the MCI routine 
are given in Tables 1-3. The column densities in Table 2 
are calculated with the upper limits of the covering fac- 
tors, i.e. real column densities may be up to 50% higher if, 
for instance, C ~ 0.6. The gas density of 5 x lO^"* cm^'' , 
and the linear size of 25 pc given in Table 1 are estimated 
with the mean intergalactic value of Jgi2- However, for an 
absorber located close to the quasar the intensity of the 
ionizing continuum Jgi2 will be enhanced making the cor- 
responding gas density higher and the linear size smaller. 
For example, J912 in the vicinity of the associated system at 
0abs = 2.1470 towards HE0141-3932 (Sect. 2.2.1 in Paper I) 
was enhanced by two orders of magnitude as compared to 
the mean intergalactic level. With such an intensity the lin- 
ear size of the present absorber would be of a sub-parsec 
scale. 

In spite of all uncertainties inherent to the present sys- 
tem, the relative abundance ratios retain a quite stable pat- 
tern (Table 3). Namely, the abundances of carbon, magne- 
sium and aluminium are strongly oversolar reaching 2.5 — 3 

solar values whereas silicon and nitrogen are considerably 
underabundant compared to carbon (by 0.2 — 0.4 dex de- 
pending on the shape of the ionizing spectrum and the 
adopted covering factors for silicon and nitrogen ions). 

The derived underabundance of nitrogen to carbon, 
[N/C] < —0.2, at highly oversolar metallicity does not 
comply with models of chemical evolution of gas in the 
central parts of QSO-hosting galaxies (Hamann & Ferland 
1999) and models of galactic chemical enrichment (Calura 
& Matteucci 2004) all of which predict significant overabun- 
dance of nitrogen relative to other elements at metallicities 
above solar. However, nitrogen measurements in planetary 
nebulae show that even at high metallicities nitrogen can 
be both under- and overabundant (AUer & Czyzak 1983; 
Perinotto 1991; Richer & McCall 2008). Low ratios [N/C] 
< at near-solar metallicities were reported also for inter- 
galactic absorbers by Reimers et al. (2005) and Jenkins et 
al. (2005). 
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Table 2. Calculated column densities for atoms and ions identified in metal-rich absorbers (for references, see Table 1). 



Redshift, 


Hi 


Cii 


Nil 


Oi 


Mg II 


Sin 




Cm 


Niii 


Ovi 


Fell 


Sim 






Civ 


Nv 




Aim 


Si IV 



2.5745* 
1.7529* 
1.5080* 



2.4158* 
A 



B 
C 

1.7315 

1.6514 
1.5770 



2.8980* 
1.7963 

2.3520 



(9.5±0.3)E14 
(1.0 - 1.5)E16 
(1.0 - 1.5)E16 

(2.2±0.1)E13 
(7.5±0.5)E15 
(1.0±0.2)E13 
(0.9 - 1.2)E16 

(3.5 - 4.5)E15 
(0.9 - 1.0)E14 



1.7817 (2.5 - 3.0)E15 



(1.1 - 2.3)E13 
(0.9 - 1.2)E16 

(0.7 - 1.3)E13 



HEI347-2457 (z^m = 2.578) 



(1.6±0.2)E13 
(1.9±0.3)E14 
(3.2±0.1)E13 
(1.5±0.1)E14 

(3.2±0.5)E14 



(1.8±0.2)E13 
<5.0E11 



<1.0E13 



(1.7±0.2)E14 <7.0E12 

HEOl 51-4326 (zem = 2.775) 



<1.1E13 
< 1.5E12 
<4.0E12 
<7.0E12 



< 5.0E11 
(1.5±0.2)E13 
(1.2±0.1)E13 

< 1.0E12 
(6.5±1.0)E12 

<1.0E12 — 

(6.5±0.5)E12 *(5 - 6)E12 < 4.5E12 

(7.0±0.5)E12 — 

(1.0±0.1)E14 — — 

(3.2±0.5)E14 <1.5E13 

HEOOOl-2340 (Zen, = 2.26) 

*3.7E13 — <1.0E12 

(4.7±0.2)E13 — 

(7.6±0.7)E12 — — 

(5.2±0.5)E13 t(0.8 - 1.1)E13 

HEOI4I-3932 ("Zem = 1.80) 

(1.64±0.08)E13 — — 

(2.2±0.2)E13 <1.4E12 

HE2347-4342 ('Zem =2.88) 

<3.0E11 — — 
(5 - 9)E12 



(1.2±0.1)E14 
t5.9E13 

(1.2±0.1)E13 



(2.6±0.2)E13 



(2.8±0.2)E14 
(1.0±0.2)E13 



<5.0E11 
t3.6E12 
(3.2±0.1)E13 



Q0329-385 (z^^ = 2.435) 



+ 2.9E12 
(4.0±0.2)E13 



(1.2±0.1)E14 



(1.3±0.2)E12 

(3.6±0.5)E11 
(6.8±0.2)E12 

< 3.5E11 
(1.7±0.2)E12 
(1.1±0.1)E13 
(2.1±0.2)E12 
(2.2±0.2)E12 



(9.0±1.0)E12 
(3.4±0.3)E11 
(1.3±0.2)E12 



4.1±0.4)E12 
2.0±0.2 Ell 
<1.3E11 
<3.0E11 

<1.3E11 



(1.2±0.2)E12 
(1.4±0.4)E11 
(3.6±1.0)E11 



(0.9 - 1.6)E13 
'9.4±0.5)E11 
7.4±0.7)E11 



(1.4±0.2)E12 
(4.7±0.4)E12 
(3.7±0.3)E13 
(9.2±1.0)E12 
(5.5±0.5lE13 
(4.5±0.5)E13 
(1.5±0.1)E13 

(5.4±0.5)E13 



<8.5E10 

< 2.5E11 

<2.8E11 

< 2.5E11 



<2.5E11 
(1.0±0.1)E13 
(3.6±0.5)E13 
(3.8±0.5)E13 



(1.8±0.2)E12 
(2.5±0.4)E12 
(2.2±0.2)E12 

5.4±1.2)E11 
1.0±0.2)E12 

(2.8±0.3)E12 
(4.5±0.4)E12 



(1.4±0.1)E13 
(4.9±1.0)E12 
(2.9±0.3lE12 



* Incomplete covering of the background light source. 

-40 < w < 30 km s^^ , B: 130 < w < 330 km s"^ , C: 300 < t; < 370 km s'^ , see Fig. 8. 
^Deconvolved from Ly— a forest absorption. 

^Calculated with the recovered velocity and density distributions. 



Quite peculiar in the abundance pattern of the system 
under study is its low content of silicon. In optically-thin 
quasar absorbers the standard finding is [Si/C] > which 
is attributed to the SNe Type II stars as the main source 
of silicon enrichment (Songaila & Cowie 1996; Simcoe et 
al. 2006). In the Zabs = 2.5745 system the contents of C. 
Mg and Al follow the solar pattern, thus we can con- 
clude that it is not carbon which is enhanced but sili- 
con which is in fact deficient. Taking into account high 
metallicity of the absorbing gas, it is conceivable to as- 
sume that this deficit is caused by depletion of silicon into 
dust. It is well known that convective envelopes of post- 
main sequence stars (Red Giants, Red Supergiants and 



Asymptotic Giant Branch stars) as well as planetary neb- 
ulae forming around these stars are places of the effective 
dust condensation. In the oxygen-rich environment, C/0 
< 1, the main dust species are iron and silicates consist- 
ing of olivine Mg2a;Fe2(i_j:)Si04, pyroxen Mga^Fei-^jSiOs 
and quartz Si02 (Gail & Sedlmayr 1999; Jeong et al. 2003; 
Ferrarotti & Gail 2006). Proportions of olivine, pyroxen 
and quartz among the silicate dust grains are unknown, 
thus magnesium could also be depleted along with silicon 
and oxygen. However, there are reasons to suppose that 
depletion of magnesium — if any — is in fact much less 
pronounced than that of silicon (see Sect. 4.2). 
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Table 3. Abundance patters for the absorbers hsted in Table 1. [X/Y] means \og{Nx/NY) — log(iVx/iVY)©. Point estimates have 
uncertainties of 0.05 dex. 



# 


^abs 


[C/H] 


[N/C] 


[O/C] 


[Mg/C] 


[Al/C] 


[Si/C] 


[Fe/C] 


1 


2.5745 


0.4 - 0.5 


-0.2 




- 0.1 


0.0 


-0.3 




2 


1.7529 


0.1 - 0.2 


< 0.2 




-0.1 - 0.1 


0.0 


-0.3 


< -0.25 


3 


1.5080 


0.1 - 0.2 


<0.3 




k, 0.2 


k, 0.15 


0.0 


^ 0.3 


4A 


2.4158 


0.2 - 0.5 


^ 0.5 








< -0.6 




4S 


2.4158 


-1.5 --2.5 


< 0.5 








< -0.5 




4C 


2.4158 


0.2 - 0.3 


<0.5 








< -0.6 




5 


1.7315 


0.2 - 0.35 


<0.1 




0.25 


0.0 


-0.25 


-0.2 


6 


1.6514 


0.40 - 0.45 




<0.1 


0.1 -0.15 


< -0.6 


-0.8 


^ -0.3 


7 


1.5770 


~ 1.0 


0.3 - 0.5 




< 


< 


-1.0 - -0.7 




8 


1.7817 


0.15 


<0.0 




0.15 


0.15 


-0.05 


0.15 - 0.25 


9 


2.8980 


~ 1.0 


-0.3 


0.0 - 0.5 










10 


1.7963 


0.5 




0.0 


0.2 - 0.4 


0.12 


0.0 


-0.55 


11 


2.3520 


0.3 - 0.5 


0.4 


0.0 
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Fig. 3. Same as Fig. 1 but for the z^bs = 1.7529 system (HE1347-2457). The vertical axis is normalized intensities. The zero radial 
velocity is fixed at ^: = 1.7529. 



Taking into account high metalhcity, relative element 
ratios, small linear size and proximity to the quasar, we 
can suppose that the absorber is formed by fragment (s) of 
planetary ncbula(e) or AGB-star cnvclopc(s) ejected into 
the intcrgalactic space by the quasar wind. Since gas and 
dust particles have different velocities in the stream, dust- 
gas separation occurs, and the gas carried away into the 
halo of the quasar host galaxy or even into the intcrgalactic 
space becomes dust-free. However, elements condensed into 
the dust remain deficient even in the dust-free absorber. 

3.1.2. System at Za.bs = 1.7529 

The absorber at Zabs = 1.7529 exhibits a saturated Hi 
A1215 line and strong lines of many carbon and silicon ions, 
all having very simple profile shapes (Fig. 3). The blue and 
red wings of the hydrogen line look different, and the metal 
lines are shifted relative to its center. Additionally, the flat 
bottom lines of the C IV AA1548, 1550 doublet and Sim 
A1206 show non-zero residual intensity at the line centers 



which means that the background radiation is leaking. The 
covering factor for these lines can be easily set as C(C iv) = 
0.99, C(Si III) = 0.98. Clear profiles are also available for the 
doublets Si IV AA1393, 1402 and Mgii AA2796, 2803, so it is 
possible to estimate the accurate covering factors for these 
lines as well: C(Siiv) = 0.98, and C(Mgii) = 0.98. Assuming 
that ions Cii and Sin trace the same gas as Mgii, we set 
C(Cii) = C(Siii) = 0.98. 

The incomplete covering arises probably due to micro- 
lensing of the background quasar by some intervening 

galaxy. This naturally supposes a very small size of the 
absorber. In spite of the presence of only one saturated line 
of neutral hydrogen, it is possible to estimate its column 
density with a relatively high accuracy (up to 30%) since 
the velocity dispersion of the gas can be well restricted by 
the numerous metal lines. 

Given the ions of different ionization stages of the same 
element (Cii, Civ, Sin, Sim, Siiv), the spectral shape of 
the ionizing background radiation can be restored rather 
accurately (see Fig. 4). In order to describe the observed 
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Fig. 4. Same as Fig. 2 but for the Zabs = 1.7529 system 
(HE1347-2457). 



intensities, 7a, of all ions at the same value of Uq, the ion- 
izing spectrum should have a soft slope with the index 
a = 1.7—1.8 between 1 and 4 Ryd, and a drop in the 
intensity at 4 Ryd by 3 — 4 times. 

The expected position of the neutral oxygen line Oi 
A1302.16 is blended with a strong forest absorption which 
prevents setting a limit on oxygen abundance. This un- 
known abundance represents an additional degree of free- 
dom in the procedure of the spectral shape recovering: with 
scaled down (relative to carbon) oxygen abundance the 
above conditions are fulfilled for the spectral shapes which 
are softer at E > 4 Ryd. An acceptable range for the ion- 
izing spectra is shown in Fig. 4. In all cases the carbon 
content comes slightly oversolar (1.1 — 1.5 solar values). 
The physical parameters of the absorbing gas are given in 
Table 1, column densities and the element abundances — 
in Tables 2 and 3, respectively. 

The abundance pattern is almost identical to that found 
in the previous Zahs = 2.5745 system: silicon is significantly 
depleted relative to carbon, magnesium and aluminium. 
The absolute abundance of magnesium depends on the 
depth of the break at 4 Ryd: more softer spectra give lower 
Mg content, so that the ratio [Mg/C] can be by ~ 0.1 dex 
both below and above zero. A clear continuum window at 
the expected position of the Fe ii A2600 line allows us to set 
an upper limit on the iron content: [Fe/C] < —0.25, i.e. iron 
is clearly underabundant relative to carbon. Unfortunately, 
both nitrogen lines Nv AA1238. 82, 1242.80 are contami- 
nated by forest absorption, and only a not very instructive 
upper limit on the nitrogen content can be set, [N/C] < 0.2. 

The high metallicity and the relative ratios of C, Si, Mg, 
Al, with silicon and probably magnesium depleted com- 
pared to carbon and aluminium, suggest again (as it was 
for the previous Zabs = 2.5745 absorber) the dust-forming 
envelopes of the post-main sequence stars or planetary neb- 
ulae as a possible source of the observed absorption. The 
underabundance of iron relative to carbon can be explained 
by different reasons: it may be intrinsic due to a relatively 
low contribution from SNe la to the enrichment of gas, or it 




Relative Velocity (l^rn ; 

Fig. 5. Mgii in the Zabs = 1.7529 system (HE1347-2457). Panel 
a: observed (points) and synthetic (solid lines) profiles when 
Mg II lines arc centered using the solar composition of the iso- 
topes 2*Mg:2^Mg:2''Mg = 79:10:11. Panel b: observed (points) 
and synthetic (solid lines) profiles when Mg II lines are shifted by 
0.6 km s~^ which corresponds to the inverse isotope composition 
^*Mg:^^Mg:^^Mg = 11:10:79. The corresponding signal-to-noise 
(S /N) mean value per resolution clement at the continuum level 
and the obtained Xmin per degree of freedom are shown. 



may be produced in the AGB-stars themselves through the 
depletion onto dust grains and/or s-process (Hcrwig 2005). 

Clear and strong lines of Mg ii make it possible to con- 
duct one more test whether the AGB-stars are at play or 
not. Theory predicts that in the AGB-stars with masses 
M > AMq the heavy magnesium isotopes ^^Mg and ^^Mg 
should be enhanced (Karakas & Lattanzio 2003). Solar iso- 
topic composition is 24Mg:25Mg:26Mg = 78.99:10.00:11.01, 
but observations of red giant stars in globular clusters show 
in some cases much higher input of ^'''Mg and ^^Mg (Yong 
et al. 2003, 2006). Due to secondary character of ^^Mg and 
^^Mg the content of these isotopes is expected to scale with 
metallicity. 

The wavelengths of the resonance transitions of Mgii 
are 2796.3553 A for ^'^Mg, 2796.3511 A for ^SMg, and 
2796.3473 A for ^SMg (Morton 2003). For solar isotopic 
composition Morton gives the weighted mean value Arest = 
2796.3543 A. If the abundances of the heavy isotopes are en- 
hanced, then the wavelength Arcst becomes smaller, i.e. the 
observed at a given redshift Mg ii lines are shifted toward 
lower radial velocities as compared with their positions ex- 
pected for the solar isotopic composition ratio. 

In our approach to calculate synthetic profiles we as- 
sume that all lines are produced by the same absorbing gas 
with fluctuating density and velocity fields. Profiles of spe- 
cific ions may look different due to different responses of 
ions (ion fractions) to the local gas density which in the 
optically thin case are determined entirely by the spectral 
energy distribution of the ionizing background. This means 
that profiles — in spite of being difFcrcnt should retain 
definite consistency governed by both the distributions of 
gas density and velocity along the line-of-sight and by the 
ionizing background. However, with magnesium lines cen- 
tered using the solar composition wavelengths the synthetic 
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Fig. 6. Same as Fig. 1 but for the ^abs = 1.5080 system (HE1347-2457). 



profiles of Mgii and other low- ionization transitions Cii 
A1334 and Sin A1260 look a bit incoherent (Fig. 5a). 

On the contrary, when the observed profiles of Mg ii are 
shifted by, e.g., 0.6 km s^^ towards higher radial velocities 
(which would correspond to the inversion of the solar ra- 
tio of 24Mg to 26Mg, i.e. 24Mg:25Mg:26Mg = 11:10:79), the 
quality of the fitting improves considerably (Fig. 5b). 

However, with the available spectral data this test is 
to a certain extent qualitative: positive shift of magnesium 
lines is unambiguously preferred, but the accurate value of 
this shift cannot be determined since the fittings with Mg II 
lines shifted by 0.5 km s~^ and by 0.7 km s~^ are statis- 
tically indistinguishable. In fact, accurate measurements of 
the isotopic composition by means of a line shift require 
a specially observed and processed quasar spectrum (e.g., 
Levshakov et al. 2007; Molaro et al. 2008) and cannot be 
carried out with the present spectrum of HE1347-2457. 

As already mentioned above, the linear size of the Zabs — 
1.7529 system is probably very small. The upper limit (as- 
suming the absorber is intergalactic) is L < 500 pc. The 
continuum window at the expected position of the Cii* 
A1335 line allows us to set an upper limit on its column 
density, N{Cn*) < 9.0 x 10^^ cm^^ which results in an 
upper limit on the gas density of Uh ^ 0.14 cm~^ (see Sect. 
2.2.1 in Paper I for details of calculations). With the total 
hydrogen column density of (5 — 10) x 10^^ cm~^ this gives 
L > 12 pc. 

3.1.3. System at Zabs = 1-5080 

This system exhibits one saturated and very noisy hydro- 
gen line (lies at the edge of observational wavelength range) 
accompanied by strong lines of metal ions in different ion- 
ization stages (Fig. 6). Clear central parts of the doublets 
Civ A1548,1550 and Si IV A1393, 1402 allow us to calcu- 
late accurately the corresponding covering factors: C(Civ) 
= C(Siiv) = 0.99. The same covering factors were adopted 
for all other lines. Given only one pair of the subsequent 
ions. Si II and Si iv (C ii A1334 is blended with a deep forest 
absorption and cannot be deconvolved) and unknown rela- 



tive element abundances, it is impossible to conclude about 
the SED of the ionizing radiation. However, the strong C IV 

and the pronounced Fe II lines point to a spectrum at least 
as hard at i? > 4 Ryd as that estimated from the previous 
Zabs = 1.7529 system (Fig. 4) — otherwise an overabun- 
dance of iron to carbon significantly exceeds solar value, 
[Fe/C] >1.0. 

For the same reason spectra with the indices a ~ 
1.5 — 1.8 between 1 and 4 Ryd are preferable. With such 
spectra, abundances of carbon and silicon come out to be 
slightly oversolar, however without deficit of silicon (cf. with 
column densities in the Zabs = 1-7529 system): [C/H] ^ 
[Si/H] ~ 0.1 — 0.2, whereas magnesium, aluminium and iron 
are strongly enhanced relative to carbon: [Mg/C] > 0.2, 
[Al/C] > 0.15, [Fe/C] > 0.3, with very probable overabun- 
dance of iron to magnesium as well, [Fe/Mg] > 0.1. There 
is a continuum window at the expected position of the N v 
A 1238 line which provides a conservative upper limit on 
the nitrogen abundance [N/C] < 0.3. The physical parame- 
ters, column densities and derived abundances are given in 
Tables 1, 2 and 3. Very similar abundances were previously 
obtained for the narrow-line associated system at Zabs = 
1.7817 towards HE0141-3932 (# 8 in Tables 2 and 3, see 
also Sect. 4.2 in Reimers et al. 2005). Such pattern — with 
iron enhanced relative to a-clcmcnts is often observed in 
giants in local dwarf galaxies (Venn et al. 2004; Bonifacio et 
al. 2004; Tautvaisiene et al. 2007) and was also registered 
in some metal-rich giant stars in the Milky Way (Fabbian 
et al. 2005). 

As it was the case for the 2:abs = 1-7529 system (Sect. 
3.1.2), simultaneous fitting of all lines in the present system 
results in systematically shifted low-ionization lines when 
Mgii is centered using the solar isotopic ratio (Fig. 7a), 
whereas Mgii lines shifted by ~0.6 km s^^ redward pro- 
duce almost perfect fitting (Fig. 7b). Thus, the enhanced 
content of heavy Mg isotopes seems very probable which 
means that AGB-stars contributed to the metal enrichment 
of gas in the Zabs = 1-5080 system. 
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3.2. Quasar HE0151-4326 
3.2.1. System at ^;abs= 2.4158 

The absorbing complex at Zabs = 2.4158 was firstly men- 
tioned in Aracil et al. (2004). It is spread over 350 km s~^ 
and can be divided into three parts (Fig. 8): two systems 
with weak H i Ly-a lines and pronounced metal lines of 
C III A977 and C iv AA1548, 1550 centered at z; = km s'^ 
(subsystem A) and v = 347 km (subsystem C, ^abs = 
2.4196), and a subsystem B comprising a strong H I absorp- 
tion in the range 140 < t; < 340 km s~^ (zahs = 2.4180) 
accompanied by only one prominent metal line - C III A977. 
That the observed absorptions are due to C ill is quite cer- 
tain: lines are narrow (i.e. cannot be hydrogen lines from 
the Ly-a forest), aligned in velocity with the hydrogen ab- 
sorption, and there are no any plausible metal contaminants 
from the identified intervening systems. 

In the subsystems B and C, there are weak and narrow 
hues at the expected positions of the Niii A989 (in C - 
blended with broad and shallow Hi forest absorption). No 
candidate for blending from other systems was found, so 
quite probably these weak absorptions arc indeed due to 
Niii A989. From the Nv doublet, only a small unblended 
part of the Nv A1238 line is present in the subsystem A. 
The O VI doublet has an unblended part in the O vi A1037 
line seen in the subsystem C. 

The C II A 1334 lines coincide with a broad and shallow 
intrinsic quasar absorption of Hi Ly-a (HE0151-4326 is a 
mini-BAL quasar), but the absence of any pronounced fea- 
tures at the expected positions of the C ii lines in all three 
subsystems leads to the conclusion that these lines should 
be very weak. Continuum windows at the positions of the 
silicon lines Sim A1206 and Siiv A1393 give upper limits 
on their column densities. Unfortunately, both lines of Mgll 
A2796, 2803 coincide with strong telluric absorptions. 

With ions available in all three sybsystems it is not pos- 
sible to restore the spectral shape of the underlying ionizing 
background uniquely: the observed line profiles can be de- 
scribed with UVB spectra ranging from pure power laws 
to power laws with intensity break of different depths at 4 
Ryd. The blue wing of the H i Ly-a line in the subsystem A 
turned out to be inconsistent with the assumption of con- 
stant metallicity throughout the absorber. The blue wing 
of the synthetic Ly-a profile was then calculated on base of 
density and velocity distributions restored from metal ions 
and the red part of the H i Ly-a assuming a constant metal 
content of the absorbing gas. 

The column densities of all ions are listed in Table 2. 
The subsystems A and C reveal almost identical ratios 
Civ/Ciii ^ 1 and Civ/Hi ^ 0.5 indicating that these 
absorbers have very similar ionization states and metallic- 
ities. In contrast, the absorber B with its ratio Civ/Ciii 
< 0.15 and more than two orders of magnitude higher col- 
umn density of neutral hydrogen, iV(H i) = 7.5 x 10^^ cm~^ , 
seems to be different both in physical state of the gas, and 
in metal content. 

In spite of the uncertainty in the UVB shape, element 
abimdances retain a remarkably stable pattern (Table 3). 
In particular, in the subsystems A and C at the value of the 
ionization parameter U corresponding to Tciv/f Cm ^ 1 
(Tj is a fraction of ion i) all ionizing spectra give oversolar 
abundance of carbon which can reach a few solar values 
(the harder at E > A Ryd spectrum the higher carbon 
content). On the other hand, the fractions of silicon ions 
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Fig. 7. Same as Fig. 5 but for the Zabs = 1.5080 system 
(HE1347-2457). In contrast to Fig. 5a where both lines Sin and 
Mgii are equivalently strong, here Fell is significantly weaker 
than Mgii and, hence, the weight of Mgii dominates in the 
minimization of the objective function leading to a prominent 
shift of the Fe ii line. 

Sim and SiiV at U corresponding to Tciv/Tcm ~ 1 are 
such high that only a significantly undcrabundant silicon 
can comply with upper limits on A^(Silll) and A^(Silv): 
[Si/C] < —0.6. At this U for all spectra we also have Tc m = 
Tniii- If the line in the subsystem C is indeed due to Nlll 
A989 absorption, then we obtain a clear overabundance of 
nitrogen to carbon, [N/C] ~ 0.5. 

Unfortimatcly, the upper limit on N{0 vi) from the sub- 
system C cannot be translated into a reasonable boimd on 
the oxygen abundance: the O VI fraction is extremely sen- 
sitive to the energy level above 4 Ryd, but since this level 
cannot be estimated from the available ions, variations of 
the oxygen abundance for applicable ionizing spectra ex- 
ceed an order of magnitude. 

In the subsystem B, the available limits Civ/Ciii < 
0.15 and Nll/Nlli < 0.6 provide bounds to the acceptable 
U range: -3.2 < log(C/) < -2.25. In this U range, log(THi) ~ 
-2.2 - -3.2, the fractions of Cm and Niii remain almost 
constant, log(Tciii) ^ log(TNiii) ~ -0.05, and the fraction 
of Sim changes from -0.1 to -0.2. All values depend only 
weakly on the spectral shape. Thus, the measured column 
densities lead to the following abundances: [C/H] ~ -2.5 - 
-1.5, [Si/C] < -0.5, and [N/C] - 0.5 - 0.6 (providing that 
Nlll A989 absorption in subsystem B is real). 

This kind of the abundance pattern when a low carbon 
content (1/300 to 1/30 solar) is accompanied by a strong 
underabundancc of silicon (and probably by 3-times over- 
abundance of nitrogen) is highly unusual for a stand-alone 
intergalactic absorber with a comparable neutral hydrogen 
column density, A'^(H i) ~ 10^^ cm~^ , but it coincides with 
[Si/C] (and [N/C]) measured in the subsystems A and C. 
Such a similarity points to a physical connection between 
all three absorbers. Note that the ionization parameter Uq 
in the metal-poor subsystem B is signifucantly lower than 
that in the metal-rich subsystems A and C, i.e. the gas den- 
sity in the subsystem B is higher (the kinetic temperature 
is higher as well because of low metallicity). 
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Fig. 8. Same as Fig. 1 but for the complex absorption system at «abs = 2.4158 towards HE0151-4326. Three subcomponents cover 
the ranges: -30 < ?; < 40 km s"^ (A), 140 < u < 340 km s"^ (B), and 300 < u < 380 km s"^ (C). 



Absorption systems with strong metal lines and low col- 
umn densities of neutral hydrogen, A/'(Hl) < 10^^ cm~^ , 
such as the subsystems A and C, arc known from the liter- 
ature (e.g., Schaye et al. 2007). For these systems the de- 
rived metal abundances may be; artificially boosted because 
of a possible overionization of Hi (see, e.g., Sect. 2.3.1 in 
Paper I). However, the upper limit on A''(0 vi) < 4.5 x 10^^ 
cm ~^ and the low value of the ratio iV(Ovi)/A''(Civ) < 
0.65 in the subsystem C show that gas is close to the ion- 
ization equilibrium. Thus, the measured high metal abun- 
dances are probably real. 

The obtained abundance pattern with silicon deficient 
(and nitrogen overabundant) to carbon, resembles patterns 
(taking silicon as a proxy to oxygen) observed in plane- 
tary nebulae in the Milky Way as well as in the LMC and 
SMC (StanghcUini 2007). Carbon (and nitrogen) arc obvi- 
ously enhanced due to dredge-up processes in the progeni- 
tor star(s). Whether a low content of silicon is intrinsic or 
affected by some depletion into dust is not clear since there 
are no other elements available. 

It is known that many planetary nebulae/AGB-stars 

have cold outer envelopes (sec, e.g., Marengo 2009, and ref- 
erences therein). In order to explain systematically higher 



abundances derived from weak optical recombination lines 
as compared to the abundances resulting from collision- 
ally excited lines, a model of planetary nebula (PN) as- 
suming the presence of H-deficicnt (i.e. metal-rich) inclu- 
sions embedded in the diffuse nebular gas is discussed (e.g., 
Middlemass 1990; Liu et al. 2004; Wesson & Liu 2004; 
Zhang et al. 2005). Accounting for this model, the ^^abs = 
2.4158 absorption system can be interpreted as PN frag- 
ment (s) transported into the IGM by the AGN/galactic 
wind. 

In Paper I, we described a system at 2;abs = 2.3520 to- 
wards Q0329~385 with an oversolar carbon content and 
overabundance of nitrogen to carbon (Table 3). This sys- 
tem exhibits a weak H I Ly-a accompanied by strong metal 
lines (Civ, Nv, Ovi), i.e. resembles the subsystems A and 
C, but is much higher ionized. The Zahs = 2.3520 system 
is shifted by ^^7400 km s~^ from Q0329 385 and was clas- 
sified in Paper I as a probable eject from the quasar host 
galaxy. While analysing this system, we have not consid- 
ered the nearby absorption systems. However, detached by 
only 230 km s~^ (^abs = 2.3545) there is an absorber with 
properties similar to those of the subsystem B: strong Hi 
Ly-a along with clear continuum at the positions of the 
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Fig. 9. Same as Fig. 1 but for the complex absorption system at 2:abs = 2.3520 towards Q0329-385. The subsystem at v = km s 
is described in Paper I, Sect. 2.4. The subsystem at v = 230 km has ^abs = 2.3545. 
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Fig. 10. Same as Fig. 1 but for the Za.hB = 1.7315 system (HE0151-4326). 



Cii, Civ, Sin, Siiv, and Nv lines (Fig. 9). Weak absorp- 
tions are present at the positions of Cm A977 and Ovi 
A1032, but since both fall in the Ly-a forest, the observed 
intensities are only upper limits for these lines. The column 
density of H i can be estimated quite accurately (Ly-/3 and 
Ly-7 available): iV(Hi) = (4.5 ± 0.3) x lO^^ cm^^ . With 
such a high column density Ly-a should be fully saturated - 
but it shows a flat bottom with the residual intensity of 
0.025 (see insert in Fig. 9). This means that the absorber is 
small and does not completely cover the light source. At the 
same time, the system at z^hs = 2.3520 with A^(Hl) ~ lO-"^^ 
cm~^ , the total hydrogen column A'^(H) of a few units of 
10^'^ cm~^ and the gas density n > 10~* cm~^ (Table 1) 
does not show incomplete coverage, i.e. in any case it is 
not smaller than the system at Zabs = 2.3545. Therefore, 
the absorber with TV (Hi) = 4.5 x 10^^ cm~^ should have 



A''(H) < 10^* cm~^ and n > 10~^ cm~^ which is realized 
at \ogU < —2, i.e. it cannot be a highly ionized system 
(and apparent absorption at the position of Ovi A1032 is 

obviously due to some forest line). With the upper limit on 
Af(Ciii) < 3.6 X 10^^ cm~^ we obtain for the carbon con- 
tent the estimate [C/H] < —1.0. Unfortunately, positions of 
Sim A1206 and Nlll A989 arc blended and any conclusive 
estimates on silicon and nitrogen contents are impossible. 
Thus, the system at Zabs = 2.3545 comprises low-ionization 
and metal-poor gas — as was the case for the subsystem 
B, whereas the Zabs = 2.3520 absorber — just as the sub- 
systems A and C — is formed by metal-rich and highly 
ionized gas. Surprisingly, even the velocity offset between 
the Zabs = 2.3520 and Zabs = 2.3545 systems exactly corre- 
sponds to that between the subsystems A and B. 
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Fig. 11. Same as Fig. 2 but for the z^bs = 1.7315 system 
(HE0151-4326). 



3.2.2. System at -Zabs= 1.7315 

This system is very similar to that at 2;abs = 1-7529 towards 
HE1347-2457 (Sect. 3.1.2). A narrow saturated Hi A1215 

line is blended in the blue and red wings by hydrogen lines 
from the Ly-a forest (Fig. 10). Both the Siiv hues 1393 
A and 1402 A coincide with broad intrinsic Ovi absorp- 
tion (HE0151-4326 is a mini-BAL quasar), but the profile 
of Si IV A1393 can be deconvolved. The line Mgii A2796 
is contaminated by telluric absorption. The lines of C iv 
AA1548, 1550 faU in the Ly-a forest, Civ A1550 is blended 
in the blue wing. Thus it is not clear whether the incom- 
plete coverage is realized or not. In the following all lines are 
considered as having covering factors of unity. Acceptable 
shapes of the ionizing spectra arc shown in Fig. 11, the 
derived physical parameters, column densities and the cor- 
responding element abundances are presented in Tables 1, 
2, and 3. 

Silicon is considerably depleted relative to carbon, 
whereas the magnesium to carbon ratio is almost two times 
higher than the solar value. There is also deficit of iron 
with respect to carbon. Similar pattern characterized by 
the strong overabundance of magnesium and deficit of iron 
(deficit of silicon was not reliably detected), was found in 
the absorption system at 2;abs= 1-7963 towards HE2347- 
4342 (# 10 in Table 1, see also Sect. 3.1.2 in Agafonova 
et al. 2007). The overabundance of magnesium is observed 
in AGB-stars and is related to the so-called 'hot-bottom' 
burning (Herwig 2005). 

Again, as in the case of the systems at Zabs = 1.7529 and 
z^y,^ = 1.5080 towards HE1347-2457 (Sect. 3.1.2 and 3.1.3, 
respectively), the shift of the magnesium lines to ~ 0.6 
km s~^ leads to a noticeable improvement of the fitting. 
The present system shows a narrow but blended line at the 
position of Nv A1238. Since the expected position of Nv 
A1242 is blended with a strong forest absorption, the col- 
umn density obtained from the Nv A1238 line is in fact 
an upper limit. A conservative limit on the nitrogen abun- 
dance shows that the overabundance of nitrogen to carbon 
does not exceed 0.15 dex (1.4 times). 



3.3. Quasar HEOOOl-2340 

3.3.1. System at Za,hs = 1.6514 

This system consists of a narrow hydrogen Hi A1215 line, 
saturated and blended in the blue and red wings, and many 
strong lines of metal ions in different ionization stages 
(Fig. 12). The simultaneous fitting of all lines assuming 
constant metallicity throughout the absorber makes it pos- 
sible to restore the wings of the hydrogen Ly-Q line. The 
line of Cll A1334 is blended with a weak forest absorp- 
tion and can be deconvolved. The expected positions of 
the Nv AA1238, 1242 lines are contaminated with a strong 
intervening absorption which prevents setting any reason- 
able limit on the nitrogen abundance. The shift of the Mgii 
A2796; 2803 lines by ~ 0.6 km s^^ results in the improved 
fitting of all low-ionization lines, C il A1334, Si II A1526, and 
Mgii AA2796,2803 (see Sect. 3.1.2). It cannot be excluded 
that the magnesium profiles are affected by an incomplete 
coverage, C(Mgii) ^ 0.95 — 0.98, but the current quality of 
the spectral data does not allow us to make unambiguous 
conclusion. The physical parameters are given in Table 1. 
The column densities listed in Table 2 correspond to C = 1 
for all lines. 

Very similar ratios Civ/Cii = Siiv/Sill = Sim/Si IV 
~ 1.1 require a hard ionizing spectrum at 1 < < 4 Ryd 
with a ^ 0.8 — 1 and with the intensity break &t E = 4 Ryd 
of about 1 dex (Fig. 13). The resulting element abundances 
are given in Table 3. Noticeable is an extremely strong 
depletion of silicon to carbon and magnesium, [Si/C,Mg] 
~ -0.9. The hue Fell A2382 falls in the spectral region 
with many weak telluric lines, so it is not clear whether the 
absorption at the expected position of this lines is indeed 
due to Fell (no other iron lines are available). In any case, 
even an upper limit on the iron content shows that iron is 
depleted as well, [Fe/C] < —0.3. In contrast to all systems 
described above the present one has a deficit of aluminium: 
[Al/C] < —0.6. A continuum window at the expected po- 
sition of the O I A1302 line allows us to set a very conser- 
vative limit on the oxygen abundance, [O/H] < 0.5, which 
excludes any overabundance of oxygen to carbon. The lin- 
ear size along the line-of-sight is small with an upper limit 
(the intergalactic Jgi2 intensity is assumed) of L < 30 pc. 

3.3.2. System at Zabs = 1.5770 

There are three systems — at -Zabs = 1-5770, 1.5810, and 
1.5855 — located quite closely in the radial velocity range 
(A?; ^ 1000 km s"-'^ ) and showing a similar set of strong 
and narrow lines of different metal ions. However, in the two 
latter systems most lines are blended, so the quantitative 
analysis is possible only for the former absorber which is 
shown in Fig. 14. 

From silicon lines, only Sim A1206 and Siiv A1393 are 
present, Siiv A1402 is blended with deep forest absorption. 
Both lines fall in the Ly-a forest range and strictly speak- 
ing should be considered as upper limits (although no metal 
contaminants were found) . A narrow absorption at the posi- 
tion of N v A 1238 is very probably due to this ion (no metal 
candidates for blending), its red wing is blended with shal- 
low forest line, but can be deconvolved. Mgii A2796 coin- 
sides with a telluric line, weak absorption at the position 
of Mgii A2803 can also be telluric. 
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For any type of the tried ionizing spectra the carbon 
content comes out extremely high — above 10 solar values. 
The apparent profile of hy-a is inconsistent with the as- 
sumption of constant metallicity throughout the absorber 
(see synthetic profiles in Fig. 14). Note that similar in- 
consistency between profiles of hydrogen lines and hncs 
of metals was also detected in Paper I for two absorbers 
with very high metallicities: the system at Zahs = 2.898 to- 
wards HE2347-4243, and the system at ^;abs= 2.352 to- 
wards Q0329-385. Both of them are associated. 

The physical parameters and column densities arc given 
in Tables 1-3. At such high metallicities as obtained for the 
2abs= 1.5770 system the dependence of ion fractions on 
the absolute metal abundances (which are unknown) over- 
rides the dependence on the SED, so it is not possible to 
bound the acceptable range of spectral shapes. However, 
UVB spectra which are softer at > 4 Ryd than the model 
spectrum for z ^ 2 of Haardt&: Madau (1996) should be 
probably excluded since at U corresponding to the observed 
ratio Cii/CiV they do not reproduce the observed upper 
limits on Si ill and Si IV and give an extremely large (more 
than one magnitude) overabundance of nitrogen to carbon. 
The previously described system at ^;abs = 1-6514 for which 
the spectral shape can be restored is detached in the ve- 
locity space by ^8000 km s~^ . In principle, such velocity 
dispersions could be possible if absorbers were ejected by a 
quasar located transverse to the line-of-sight. In any case, 
for all types of trial spectra silicon remains considerably 
underabundant to carbon, [Si/C] < —0.7. Also for all spec- 
tra, the deconvolved N v A1238 line leads to overabundance 
of nitrogen to carbon. The value of [N/C] depends on the 
softness at E > 4 Ryd: the softer the UVB spectrum the 
higher [N/C]. 

Extreme high abundances of carbon and nitrogen might 
be artificial because of non-equilibrium ionization (H i ove- 



rionized) in the system considered. However, more plausible 
explanation seems to be a really H-deficient gas from the 
inner parts of the envelopes of the AGB-stars. The low ra- 
tio of silicon to carbon and nitrogen, [Si/C,N] < —1, can 
be explained by the enhancement of C and N due to the 
dredge-up processes and by the depletion of silicon into 
dust. 



4. Discussion 

4.1. Abundance patterns of high metallicity gas 

In the present paper we considered 7 optically-thin systems 
with solar to oversolar metallicities. Thus, together with 
previously described absorbers, total sample of metal-rich 
systems consists of 11 absorbers with redshifts 1.5 < z < 
2.9 identified in spectra of 6 quasars. To simplify the read- 
ing, in the following the systems will be referred to by their 
numbers as given in Tables 1 and 3. Metallicity is indicated 
on base of the measured abundance of carbon. 

The majority of the systems except # 9 (Tables 2, 3) 
show abundance patterns which — in spite of being dif- 
ferent — can nevertheless be related to outflows from low 
and intermediate mass stars {IMq < M < IOMq) in their 
post-main sequence evolution stages. 

The observed diversity of patterns supports the assump- 
tion that most metal-rich absorbers have very small sizes 
{L ~ 10—500 pc) as compared with the intergalactic clouds, 
i.e. absorbers are formed by outflows from only a few stars 
(or even from a single one). 

The system # 9 shows extreme high metal content ex- 
ceeding 10 solar values. It is not clear whether the appar- 
ent metallicity is artificially boosted due to the overion- 
ized state of H I or it is intrinsic and thus indicating re- 
ally hydrogen-deficient gas (see Sect. 2.3.1 in Paper I). In 
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Fig. 13. The ionizing spectrum (solid line) corresponding to the 
ionization state observed in the Zabs ~ 1.65136 system (HEOOOl- 
2340). The intergalactic ionizing spectrum at 2 = 1.6 modeled 
by Haardt & Madau (1996) is shown by the dotted line. The 
spectra are normalized as in Fig. 2. 

any case, the overabundance of oxygen to carbon points 
to SNe II as the main source of gas enrichment. The pres- 
ence of only one SNe Il-enriched system in our sample is 
probably a selection effect: in quasar spectra there is a 
non-negligible population of absorbers revealing strong C IV 
lines along with weak Hi absorption (sec, e.g., Schaye et al. 
2007). These systems may indeed contain high metallicity 
gas enriched prevalently by SNe II, but an accurate quan- 
titative analysis is in most cases prevented by the absence 
of lines in subsequent ionization stages. 

Fig. 15 shows the ratio [N/C] as a function of the car- 
bon abundance [C/H] which in the present case serves as 
an indicator of the overall metallicity. No positive correla- 
tion of [N/C] with metallicity is seen: at any value of [C/H] 
nitrogen can be either over- or under abundant. Similar find- 
ings were reported for planetary nebulae (AUer & Czyzak 
1983; Perinotto 1991; Richer & McCall 2008) and for high- 
metallicity gas near AGN (Fields et al. 2005). This indicates 
that the nitrogen enrichment occurs quite irregular and is 
probably governed by many processes. On the other hand, 
the spread of the nitrogen contents can be considered as an 
additional argument in favor of a very small linear size of 
the metal-rich absorbers (formed by outflows from only a 
few stars). 

Absorbers # 1, 8, 9 and 11 belong to the so-called asso- 
ciated, or proximate systems, i.e. they are located close to 
the background quasar and are formed in gas ejected from 
the host galaxy. All other systems are formally interven- 
ing. However, they demonstrate same physical properties 
as the above associated — compare, for example, # 1 and 
2, # 8 and 3, # 11 and 4. Two conclusions follow. Firstly, 
many metal-rich intervening systems originate in gas ex- 
pelled from the quasar host galaxies located transverse to 
the line of sight. This is also supported by the direct ob- 
servations of HE2347-4342 where transverse quasars were 
detected just at the redshifts of such systems (Worseck et al. 
2007). Secondly, associated systems not ncssesary originate 
in gas ejected from the circumnuclear regions — they can be 



interstellar objects (fragments of star/PN envelopes, diffuse 
clouds etc.) entrained by the quasar wind and transported 
into the intergalactic space. However, metal abundances ob- 
tained from the associated systems are often used to verify 
models of quasar chemical evolution (e.g. D'Odorico et al. 
2004; Fechner & Richter 2009). Obviously, before making 
this, the nature of a particular associated system should be 
clarified. 

4.2. Depletion of silicon into dust 

In Sect. 3.1.1 we argued that the deficit of sihcon to carbon 
obtained for the metal-rich gas may be due to condensation 
of silicon into dust in the envelopes of post-main sequence 
stars and the subsequent gas-dust separation in the stel- 
lar/galactic wind. Here we consider some aspects of silicon 
depletion into dust in more detail. 

In the oxygen-rich environment, C/0 < 1, the main 
dust species are iron and silicates consisting of olivine 
Mg2xFe2(i_x)Si04, pyroxen Mg^^Fei-a^SiOs and quartz 
Si02 (Gail & Sedlmayr 1999: Jeong et al. 2003; Ferrarotti 
& Gail 2006), with olivine and pyroxen observed mostly in 
iron-free (i.e. a; = 1) forms (Jaeger et al. 1998; Waters et al. 
1998; Draine 2003; Gutenkunst et al. 2008; Perea-Calderon 
et al. 2009). Although processes of dust growth are far from 
being clear in full detail, all reactions related to the sili- 
cate dust formation require the presence of the molecule 
SiO. This molecule is indeed observed in the envelopes of 
oxygen-rich AGB- and RSG-stars (van Loon et al. 2008; 
Sloan et al. 2008). Contrary to the previous standpoint that 
SiO polimerizes at low temperature ^600 K and therefore 
cannot form so-called seed nuclei (Jeong et al. 2003), recent 
results have shown that SiO nucleates in fact at significantly 
higher temperature T ~ 1200 K and forms small SiO grains 
over which the other dust components start to condense at 
continuously lower T (Nuth & Fergusson 2006). Thus, at 
the initial stages of dust formation the most depleted ele- 
ments are expected to be oxygen and silicon, with magne- 
sium and iron following at the later stages. Observational 
evidences of strong magnesium gradients throughout some 
planetary nebulae along with more or less uniform deple- 
tion of silicon can be considered as a support for this picture 
(Pequignot & Stasinska 1980; Harrington & Marionni 1981; 
Middlemass 1988). 

We note that systems with [Si/C] < (# 2, 5, 6 in 
Table 3) do not show deficit of silicon to iron which is not 
in line with usual assumption that iron should be stronger 
depleted into dust than silicon. In principle, it is possible 
that the depletion of iron is stronger in our systems as well: 
the intrinsic iron content is unknown and it can be enhanced 
relative to other elements as was obtained, e.g., for the ab- 
sorbers # 3 and 8 (and is observed in giant stars in local 
dwarf galaxies — see Sect. 3.1.3). However, the growth of 
dust grains includes many stages and is strongly affected by 
the environmental conditions such as gas metallicity, tem- 
perature, star pulsations etc. Under some conditions iron 
becomes more depicted than other elements, under other — 
does not. Model calculations of the dust growth processes 
in AGB-stars by Ferrarotti & Gail (2006) show that the 
relative depletion of Si and Fe depends on the metallicity 
and initial stellar mass, with silicate dust prevalent over 
iron dust in high-metallicity (solar to oversolar) AGB-stars 
with masses M < 2Mq or in massive AGB-stars M > AMq 
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Fig. 14. Same as Fig. 1 but for the Zabs = 1.5770 system (HEOOOl-2340). Identification of the Mgii A2803 line is putative since 
the stronger component Mg ii A2796 is blended with telluric absorption. 



experiencing hot-bottom burning. The hot-bottom burning 
is considered also as the cause of the enhanced content of 
heavy magnesium isotopes ^^Mg and ^^Mg (Herwig 2005) 
indications of which were detected in the absorbers # 2, 
3, 5, and 6 (Sect. 3.1.2, 3.1.3, 3.2.2, 3.3.1, and 3.3.2, re- 
spectively). Thus, we may conclude that the absorbers in 
question are formed mainly from gas expelled from massive 
AGB-stars. 

Note in passing that among metal-rich DLAs and LLSs, 
there are objects with both low [Si/Fe]: [Zn] = 0.28, [Si/Fc] 
< -0.5 (Pcroux et al. 2008), [Zn] = 0.25, [Si/Fe] < -0.4 
(Meiring et al. 2007), and high [Si/Fe]: [Zn] > 0.86, [Si/Fe] 
= 0.75 (Meiring et al. 2008), [Zn] = 0.12, [Si/Fe] = 0.7 
(Prochaska et al. 2006). Additionally, many DLAs shows 
differential depletion into dust with values [Si/Fe] vary- 
ing an order of magnitude from one component to another 
(Centurion et al. 2003; Peroux et al. 2006; Quast et al. 
2008). 



4.3. Evolutional effects in distribution of metal- rich absorbers 

On base of the observed lines, the metal-rich absorbers 

can be divided into two distinct groups: systems revealing 
strong lines of both low (Cii, Mgii, Sill) and high (Cm, 
Niii, Aim, Sim, Civ, Siiv) ionization species (# 1, 2, 3, 5, 
6, 8, 10) and systems without or with weak low-ionization 
lines (# 4, 7, 9, 11). The first group has the mean ion- 
ization parameter \ogU < —2, whereas the second one is 
characterized by \ogU > —1.5. In Fig. 16, the J7- values are 
plotted against the redshifts of the systems for which they 
are obtained. The absence of systems from the second group 
at 2; < 2.3 in our sample is explained by pure selection ef- 
fects: the important for these systems lines of C ill A977, 
Niii A989, and Ovi AA1032, 1037 are shifted at z < 2.3 
to the unobservable or very noisy wavelength range, and 
the remaining Hi A1215, C iv AA1548, 1550, and eventually 
Nv AA1238, 1242, and Siiv AA1393, 1402 do not allow to 
fix the ionization parameter and, hence, to derive accurate 
element contents. 



As for low-ionization systems, we find a single system 
at z > 2 1 with Zabs = 2.5745) and six systems at z < 2. 
In contrast to the high- ionization group, this result is not 
affected either by selection criteria for absorbers or by the 
wavelength coverage of available quasar spectra (low red- 
shift boundary Zabs = 1-5 is determined by the optical limit 
to observe Hi A1215) and, thus, refiects a real distribution 
of \ow-U absorbers over redshift. 

The epoch z ~ 2 is known to represent a peak activity 
in both the cosmic star formation rate and the assembly of 
galaxies and it is also characterized by the maximum space 
density of the luminous quasars (Fan et al. 2001, Dickinson 
et al. 2003, Chapman et al. 2005, Rudnick et al. 2006). 
Galaxies at z ^ 2 reveal star formation rates of 10 — 3OOM0 
yr-i and stellar masses 10*^ - 10"-'^Mq (Genzel et 

al. 2008, and references cited therein). The estimates of 
the comoving number density of the star-forming galaxies 
from different surveys (Chen et al. 2003, Adelberger et al. 
2005, Sumiyoshi et al. 2009, Cristobal-Hornillos et al. 2009) 
provide a density ng of a few 10~^h^ Mpc~^. 

Taking into account the detected jump in the number 
of low-ionization absorbers at Zabs < 2, it is interesting to 
compare the space density of these absorbers to the space 
density of galactic population at z ~ 2. 

At z = 2, the Hi Ly-a is shifted to 3647 A. Among 
the studied spectra of 20 QSOs, only half of them cover 
the wavelength range A < 3647 A, thus allowing us to 
identify low-ionization absorbers at z < 2. In spectra with 
the coverage A < 3647 A, the probability to find at least 
one low-ionization metal-rich absorber clearly exceeds 50% 
. Assuming that the mean linear size of the absorber is 
~ 50 pc, the chance that a random sightline passes through 
a metal-rich cloud requires the comoving number density 
of the absorbers Ua ^ 1/L'^D. With a typical distance 
lag £) ~ 10 Gpc (z ^ 1.8), the observed absorber fre- 
quency is realized if the comoving density of these sys- 
tems is ria ^ 4 X 10** Mpc"'^. Here the comoving density 
of star-forming galaxies is taken as Ug w 2.4 x 10~^ Mpc~^ 
(h = 0.7) which is in line with Chen et al. (2003). Then 
there are 2 x 10^ metal-rich absorbers in the vicinity of 
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[C/H] 

Fig. 15. Relative ratio [N/C] versus metallicity index [C/H] for 
metal-rich absorbers (see Table 3). 




Redshift, z 



Fig. 16. Ionization parameter log U estimated from metal-rich 

absorbers with redshifts 1.5 < z < 2.9 (see Table 2). The dashed 
line separates absorbers with low and high (7-values. 



each star-forming galaxy at ^ < 2. This estimate is quite 
rough because of the unknown geometric parameters of the 
clouds and the uncertainties in the galactic counting statis- 
tics. However, it corresponds by an order of magnitude to 
recent results on the Milky Way halo where a population 
of more than 10* weak metal-line absorbers located in cir- 
cumgalactic environment was detected (Richter et al. 2009). 
Richter et al. select their systems on base of the presence 
of a weak O l A 1302 line which supposes the ionization pa- 
rameter U of about \ogU ~ —3 (of, our system # 10). The 
physical parameters of these Milky Way 'gas wisps' arc very 
much alike to their high-rcdshift counterparts: the hydro- 
gen column density 7V(H i) ranges from 2 x 10^® cm^^ to 
3 X 10^* cm~^ , the hydrogen volume density Tin ~ 0.02 — 0.3 
cm~^ (if clouds are located at the distance 50 kpc from the 
Galactic center), the absorber thicknesses are L w 0.2 — 130 
pc, and the metallicities are ~ 0.1 — 1 solar. 

Finally, to estimate the mass of the clouds, we assume 
that they are homogeneous spheres of radius R. Then, 

4 / n \ / R\'^ 

Ma = -TTR^n^i^m^ « 0.1/z — Mq , 

3 Vcm-'^/ \pc J 

where /i is the mean molecular weight of the gas, and mn 
is the mass of the hydrogen atom. 

For ionized gas of solar metallicity {X -.Y : Z = 0.710 : 
0.265 : 0.025 by mass), /x = 0.62. With tin ^ lO'^ cm'^ 
and ~ 25 pc, we obtain for a single cloud a typical mass 
of Ma ^ 10A'/q, i.e. low-ionization metal-rich absorbers 
can indeed be produced by only a few intermediate-mass 
stars. For a star-forming galaxy with M* ~ 2 x IO^'^Mq, 
the total mass Mtot of all metal-rich absorbers reaches 
Mtot -O.OlAf,. 

The peak star forming activity at ^ ~ 2 is usually ex- 
plained by the presence of large amounts of gas supplied 
either by mergers or by smooth gas accretion onto galaxies 
along the large-scale filaments, with probable prevalence 
of the second mechanism (Genel et al. 2008; Shapiro et 
al. 2008; Genzel et al. 2008). Our data show that galaxies 
at z < 1.8 have already undergone stellar evolution suffi- 
cient to provide high (solar to oversolar) level of metallic- 
ity for the large population of stars. The build-up of met- 
als requires many generations of stars. Estimations from 



SINFONI survey of galaxies at 2; ~ 2 show that star- 
forming discs experienced constant star formation rates 
over at least 0.5 Gyr (Daddi et al. 2007). This means more 
than ten generations for intermediate mass stars — proba- 
bly enough time to acquire high metallicity. 

The high metal content has positive effect on the star 
mass-loss rates through driving the outflow both by radia- 
tive acceleration in resonance lines of metal ions and by 
radiation pressure on the dust, the formation of which is 
in turn boosted in the metal-rich environment. Thus, at 
high metallicity the recycling of matter occurs more rapid 
leading to constantly large amount of gas present in the 
interstellar space. On the other hand, starbursts and pos- 
sible AGNs in the central regions favor the formation of 
the large-scale winds which expel gas from galaxies. At 
some redshift, more gas is expelled as returned resulting 
in decline of both star forming and AGN activities. This 
point occurs somewhere below redshift 1.5. Unfortunately, 
wavelength coverage of optical telescopes limits redshifts 
of absorption systems suitable for quantitative analysis — 
below z = 1.5 the hydrogen Hi A1215 line is not seen 
and, hence, accurate metallicity of the absorbing gas can- 
not be determined. Commissioning of the Cosmic Origins 
Spectrograph — the high-resolution UV spectrograph with 
wavelength coverage from 1100 to 3200 A and the World 
Space Observatory (WSO)-UV will change the situation al- 
lowing to probe high-metallicity systems at redshifts below 
1.5. 

5. Summary and conclusions 

We report on the analysis of 11 metal-rich absorption sys- 
tems with redshifts 1.5 < z < 2.9 identified in spectra of 
6 distant quasars. Due to a special selection of the absorp- 
tion systems and the advanced methods to solve the inverse 
spectroscopy problem the obtained abimdance patterns can 
be considered as reliable. On base of these patterns, we 
identify types of stars responsible for the enrichment of the 
absorbing gas. The main results are as follows. 

1. The majority of the described metal-rich absorbers (10 
from 11) reveal abundance patterns which — in spite 
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of being different — can be attributed to the out- 
flows from low and intermediate mass stars (LIMS) in 
their post-main sequence evolutional stages (red giant 
branch and asymptotic giant branch). One absorber 
is enriched prevalently by SNe II. A small number of 
SNe Il-enriched objects does not represent real statis- 
tics, but is a consequence of selection criteria. Diversity 
of observed patterns points to small linear sizes of the 
metal-rich absorbers, i.e. every cloud can be formed by 
outflows from only a few (or even from a single) star(s). 

2. Several absorbers show a significant deficit of silicon 
with respect to other elements. This deficit can be ex- 
plained by the depletion of Si into dust in the oxygen- 
rich envelopes of the post-main sequence stars and the 
subsequent dust-gas separation in the galactic/quasar 
wind. 

3. If fines of Mgii AA2796, 2803 arc centered using the lab- 
oratory wavelengths which correspond to the solar iso- 
topic ratio '^'^Mg-.'^^Mg-.'^H'lg = 79:10:11, then profiles of 
all low-ionization lines (Cil, Mgil, Sill, Fell) cannot be 
fitted self-consistently with a required Xmin ~ 1- The 
quality of fitting improves significantly if centering of 
the Mg II lines occurs using slightly shorter wavelengths 
which correspond to an enhanced content of heavy iso- 
topes ^^Mg and ^^Mg in the absorbing gas. Although the 
quality of the available quasar spectra does not allow to 
estimate accurate isotopic composition of magnesium, 
the enhancement of heavy Mg isotopes in the systems 
at z^bs= 1.7529, 1.7315, 1.6514, and 1.5080 seems to 
be very probable. Such an enhancement is predicted for 
AGB-stars with masses M > AMq which experience a 
so-called hot-bottom burning. 

4. The abundance of nitrogen does not show correlation 
with metallicity: at any value of [C/H] nitrogen can be 
both over- and underabundant to carbon. This means 
that the nitrogen enrichment occurs quite irregularly 
and is governed by several processes. 

5. Among metal-rich systems, a group characterized by low 
ionization parameter, logf/ < —2.3, can be singled out. 
At 2; > 2, we find only one low-C/ system (zabs = 2.5745), 
whereas at z < 2 six such systems arc detected. This 
result is not affected by any selection criteria and re- 
flects the real redhsift distribution of the low-ionization 
metal rich systems. The comparison of the number den- 
sities of metal-rich absorbers and star-forming galax- 
ies at z ~ 2 shows that there should be ^ 10^ — 10* 
such absorbers around each galaxy. This coincides with 
the number of small absorbing clouds detected recently 
around the Milky Way. The redshift 2; ~ 2 is charac- 
terized by peak values of the star-formation rates and 
the quasar luminosity function which arc currently ex- 
plained by the presence of large amounts of gas in the in- 
terstellar space of high-redshift galaxies due to mergers 
and rapid gas accretion along the large-scale filaments, 
with the second option more preferable by recent ob- 
servations of the star-forming galaxies at z ~ 2. A high 
metallicity can be another cause for the peak activities 
at this redshift: large amount of metals in the envelopes 
of LIMS boosts the mass-loss from stars both due to ra- 
diative acceleration in resonance lines of metal ions and 
by radiation pressure on the dust grains. This leads to 
a rapid recycling of matter and, hence, to a constantly 
high amount of gas present in the interstellar space. 
When the amount of gas expelled from a galaxy due to 



the large-scale wind exceeds the amount of gas supplied 
by accretion and stellar winds, the star-forming and 
AGN activities decline. In order to study these processes 
through the analysis of metal-rich absorbers which 
turned out to be quite reliable indicators of cosmic evo- 
lution — joint optical/UV observations of quasar spec- 
tra are needed. This is expected to occur with future 
Cosmic Origine Spectrograph and WSO-UV. 
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